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Using coarse-grained molecular-dynamics simulations, a generic styrene-6tocfc-butadiene-6/ocA;- 
styrene (SBS) triblock copolymer under lamellar conformation is used in order to investigate the 
mutual entanglement evolution when a structure of alternating glassy (5')/rubbery (B) layers is 
submitted to an imposed deformation. By varying the amount of loop chains between each phase, i.e. 
noncrossing chains, it is possible to generate different types of S/ B interface definitions. A specific 
boundary driven tensile strain protocol has been developed in order to mimic "real" experiments and 
measure the stress-strain curve. The same protocol is also applied to a reference state consisting 
in a directed glassy homopolymers, as well as to an isotropic glassy polymer. The evolution of 
initial mutual entanglements from the undeformed samples during the whole deformation process is 
monitored. It is shown for all considered systems that initial entanglements mostly participate to 
the preyield regime of the stress-strain curve and that this network is debonded during the strain- 
hardening regime. For triblocks with a non-null amount of crossing chains, the lower the amount is, 
the longer the memory effect of the initial entanglement network in the postyield regime is. On the 
fly distributions of entanglements, which depart from the postyield regime, depict memory effects 
and long time correlations during the strain-hardening regime. For triblocks, loop chains reinforce 
these effects. 

PACS numbers: 61.41.-|-e Polymers, elastomers, plastics; 82.20.Wt Computational modeling, simulations; 
82.35 Jk Copolymers, phase transitions, structure; 82.35.Lr Physical properties of polymers. 



I. INTRODUCTION 

The macroscopic properties of multiphase polymers in- 
volve more than the simple addition of the individual 
properties of each of their constituents. They are the re- 
sult of alchemy between the material morphology (phase 
repartitions and regularity of this repartition), the prop- 
erties of each constituent, and the characteristics of the 
interfaces or the manner continuity between phases is as- 
sured. In styrene-6/ocfc-butadiene-Wocfc-styrene [SBS) 
triblocks copolymers of glassy /rubbery/glassy type un- 
der lamellar morphology [l|, Q , chains are indeed consti- 
tuted by two phases linked by covalent bonds, with each 
phase having its own mechanical properties. The way the 
phases are "connected" at the microscopic level, i.e., at 
their interfaces, may thus have an impact on the energet- 
ically involved transfer mechanisms at larger scales. For 
instance, a chain effectively forming a bridge between two 
S domains will transfer the stress more efficiently than a 
chain forming a loop between S and _B, except if the be- 
havior of the latter is modified by the presence of trapped 
entanglements. 

Coarse-grained Molecular Dynamics (MD) simulations 
are well suited to bridge the scales between global me- 
chanical properties of the material and those of each 
phases, including coarse-grained scale details for the na- 
ture of their interfaces. At global scale, the main fea- 
tures of the mechanical behavior of lamellar SBS ma- 
terials under elongation can be summarized in a linear 
elastic response with significant rigidity at low deforma- 
tion up to a yield point, a yielding and drawing at higher 
deformations with possible alteration of lamellar mor- 
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phology into a chcvronlike one [T|, and the possibility 
to reach very large elongations. Many other features can 
also be noticed concerning cyclic strain deformations and 
unloading Large deformations have a strong impact 
on the morphology of thermoplastic elastomers, affecting 
the nanostructures. In SBS triblocks, several studies 
have pointed out the presence of spherical polystyrene 
domains in the matrix Q 0] , playing the role of fillers [H , 
as well as the role of entanglements in this large strain 
regime. In MD simulations, such domain effects are hard 
to reproduce because one has to simulate huge systems 
with several repeated lamellar units, while entanglement 
effects are well accessible. 

In giassy/rubbery/glassy triblocks, the behavior of en- 
tanglements is quite complex to relate to. In the rub- 
bery phase, these ones are transient and may be mapped 
onto slip-links that slide along the chains during defor- 
mation, then allowing chains to explore many conforma- 
tions. In contrary, in the glassy domains, entanglements 
arc trapped and chains slow down and may visit differ- 
ent conformations under strong enough applied deforma- 
tion. This leads to stress dissipation and activated disen- 
tanglements. Additionally, in realistic triblocks, one has 
to take into account the contribution of loop chains, as 
well as diblock chains, for which the dynamics and stress 
transmission under deformation may be altered by the 
presence of frozen entanglements. 

In this paper, we analyze the evolution of entangle- 
ments in SBS triblocks and polymer glasses, when these 
materials are submitted to a macroscopic deformation. 
Unlike theoretical models [25l - [27| that take as input time- 
and deformation-independent entanglement sets, we an- 
alyze how entanglements evolve and depart from this as- 
sumption. Paper is then organized as follows. In Sec. HIl 
technical details about simulation model and methods 
are first given. The primitive path analysis [S^H^ (PPA) 
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is used to identify the entanglements during the whole de- 
formation process. Details about the method are briefly 
given and how this one is extended to give a statisti- 
cal representation of the entanglements. The mechanical 
response to a uniaxial deformation of the simulated sys- 
tems is discussed in Sec. lIIII Then, the Sec.|IV]deals with 
mutual entanglements statistics and how their statistical 
representation can be related to the response behavior at 
larger scales. Initial entanglements set from the unde- 
formed state are then monitored during the whole defor- 
mation process, as well as the "dynamical" or on the fly 
contribution of entanglements activated by the deforma- 
tion. Finally, in Sec. |Vl concluding remarks are given, 
and outlooks are discussed. 



II. POLYMER MODEL AND METHODS 

In the following, MD simulations are performed using a 
well-established coarse-grained model Q. The polymer 
is treated as a chain of = beads (where a 

denotes the species for block copolymers), which we refer 
to as monomers, of mass m = I connected by a spring to 
form a linear chain. The beads interact with a classical 
n-species Lennard-Jones excluded volume interaction. 
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(1) 



For particle distances, r > rc — 2.5a, the potential is 
cut-off and shifted such that it is continuous at r^,- In- 
dices {a, 13) stand for the different types of pairs, S — S, 
B — B, and S — B for SBS triblocks, and u^^ is in 
units of u'^. Adjacent monomers along the chains are 
coupled through the well-known anharmonic finite ex- 
tensible nonlinear elastic potential (FENE) in addition 
to Eq.(P, 



(2) 



UfeneW = -0.5fci?^ln 1 - (r/Ro) 



where model parameters are identical to those given in 
Ref. [1], namely k = SOw^/a^ and i?o = l-5a, chosen so 
that unphysical bond crossings and chain breaking are 
eliminated. All quantities will be expressed in terms of 
the molecular diameter a = 1. binding energy = 1 
and characteristic time tlj = \J ma? ju^ . Under such 
conditions, the glass transition temperature for the a = 
/3 = 1.0 model is sited around ksTg 0.42u°. 

Newton's equations of motion are integrated with the 
velocity- Verlet method and using a time step St = 
O.OOQtlj. Simulation cells of size (Lx, Ly, L^) are filled 
with M chains of size N. 



A. Simulated systems 

The considered triblocks copolymers are intended to 
mimic lamellar SBS triblocks such as the one depicted 
in Fig. [U in the sense of alternating glassy/rubbery lay- 
ers. Thus, our simple model consists in polymer chains 
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Fig 1: (Color online) Snapshot of a PS/PB/PS lamellar tri- 
block before deformation is applied. The glassy phases are la- 
beled with 5", while rubbery region is labeled with B. Evanes- 
cent points are used in background to represent monomers 
from each part. A loop chain (left) and a crossing chain (right) 
from the same sample are also superimposed. 



composed of 2A^5 
length per chain N 
fixed to ksT 
ing a value of u'gg 



0.3<^ 



Nb ~ 100 monomers, with a total 
= 200. The total temperature being 
the glassy phase is simulated us- 
LOu", and the rubbery phase with 
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0.5u^. In that way, the effective temperature in 
the rubbery phase is ksTs ^ Q.QvP > ksTg. Finally, 
the incompatibility between each phase is controlled by 
setting u^g = 10~^u'^, in order to maintain the lamellar 
morphology in the strong segregation regime. 

Triblocks are prepared using an enhanced version of 
the radical-like polymerization method (lol |. which is an 
extension of Gao's work We then refer to Refs. (9l.lTol| 
for further details concerning this method. 

This one has then been applied in order to generate 
four types of triblocks according to the percentage of loop 
chains, namely 0%, 30%, 60% and 100%. These loop 
chains are defined as chains that do not cross the entire 
material from a glassy phase S to the opposite S and by 
the way that also do not cross the whole rubbery phase 
B. Then, they loop back to their starting glassy phase 
while a part of the chain lies in the rubbery region. Such 
a chain is depicted in Fig. [T] 

For each system, the number of obtained chains of de- 
sired length N = 200 is respectively M = 192, M = 188, 
M = 194 and M — 191, and systems are kept un- 
der a (P) = state after a run of 10^ MD steps. 
Finally, monomers dynamic in each phases has been 
checked by computing the incoherent intermediate scat- 
tering function Fs(k, t) = {j^ J^f^i e^P {ik.Srj{t - to))) 
where Sr{t — to) = r{t) —r(to), which depicts (not shown) 
a typical glassy behavior jl2l | in S domains, and a rapid 
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relaxation behavior in the rubbery B domain. 

In order to have a reference state, a homopolymer is 
also simulated using the 0% configuration and setting 
= u'gg = 1.0u° this for all {a,/3). A subsequent 
simulation run is performed during 10'' MD steps at 
ksT = O.Su^g and (P) = using the same anisotropic 
Nose-Hoover barostat. After this run, the density of 
chain ends was checked to correspond to that of a di- 
rected homopolymer, namely, that chain ends are located 
in z ~ and z ~ L^. In the following, this directed ho- 
mopolymer will be referred to as hP. 

Finally, an isotropic homopolymer has also been gen- 
erated by annealing the previous hP polymer at ksT = 
l.Oulg and (P) = 0.5 during 2 x 10^ MD steps, and cool- 
ing it to ksT = O.Sulg for a time of 10^ MD steps. Then 
an additional run of 10^ MDS at the same temperature 
and pressure (P) = has been performed. In this case, 
the density of chain ends was checked to be spatially ho- 
mogeneous. In the following, we will refer to this sample 
as iP polymer. 

Glassy homopolymers and glassy phases in triblocks 
are then out of equilibrium. The a relaxation time Tq, 
is defined as the time at which Fs(}ii,Ta) = 1/e, where 
ki is the first peak position in the structure factor. We 
found Ta of same order than the typical time used to per- 
form the deformation experiment. In that sense, ageing 
effects may not significantly contribute to the glassy dy- 
namics upon deformation, and are expected to not alter 
the mechanical response of the materials. 



B. Boundary driven uniaxial tensile test 

A method in which the deformation is applied first at 
boundaries and then is self-transmitted to the sample 
has been developed in order to mimic experiments jll| . 
The periodic boundary condition in the z-axis (direction 
of the applied deformation) is removed, while the other 
directions are kept periodic. Hence, chains that cross 
the periodic-image boundary in the z direction were un- 
wrapped and cut. This induces length polydispersity and 
then slightly decreases the total number of monomers per 
samples. For instance, it only alters the glassy parts of 
the triblocks. Another simulation run of 5 x 10^ MD steps 
is added to ensure all samples are still under an overall 
(P) = state. 

Grips arc defined in the two glassy regions with a thick- 
ness of 2. 6055, the upper one for z > — 2.6055 and 
the lower for z < 2.6055. Then, forces and velocities 
are set to zero and vq respectively for all monomers ly- 
ing in the grips, requesting a vq ~ for lower grip. The 
lower grip is then coupled to an elastic spring with a force 
F' = —K (r^m(t) — r^„i(0)) adjusted so that the average 
lower grip velocity is null, as requested. The upper grip is 
also coupled to an elastic spring with a force also adjusted 
to preserve the requested upper grip velocity UoGz, and 
is given by - A" (r^„(t) - r^™(0) - vote^)- In both 

expressions, remit) denotes the center of mass position 
of grips at time t. 

During this boundary driven deformation protocol, a 
Berendsen thermostat is used with a heat bath oiksT = 



0.3u°^, while an anisotropic Nosc-Hoover barostat is ap- 
plied in order to keep a lateral pressure Pxx = Pyy = 0. 
The obtained deformation thus mimics a uniaxial tensile 
test. As discussed in Ref. [llj, the value of K should be 
chosen as ten times as stiff as the initial sample, in order 
to lead to a numerically stable deformation scheme; such 
a condition is fulfilled in our numerical experiment. 
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Tab I: Young's modulus Ey, hardening modulus Gr, and post- 
yield flow stress crjio-w from stress-strain curve in Fig. [21 for 
all simulated systems. Also given, the average entanglement 
length {Ne) for undeformed samples and from the primitive 
path analysis. 



C. Primitive Path Analysis 

Several theories try to explain the deformation of poly- 
mer melts in both rubbery and glassy states. In the rub- 
bery state, polymers can explore many conformations un- 
der the constraint of cross-link points (entanglements, co- 
valent bonds.. J. Entropic-based classical theory of rub- 
ber elasticity p^l - [l5{ . which neglects entanglements, may 
be used. However, in the glassy state, chains are slow 
down and may visit conformations in a noncontinuous 
way under deformation. This leads to stress dissipation 
effects [l6l - [l8l | under deformation, noisy term currently 
not present in rubber elasticity theories. 

Statistical mechanic theories improve the lack of cross- 
link effects by including entanglements to the constitutive 
relation via some additional terms [U, [2^ , as well a s by 
developing subnetwork deformation formalisms [2^ l27l |. 
It's then supposed that entanglements depart from an 
initial distribution, while nothing more is said about its 
evolution upon deformation. This is a crucial point, as 
nothing compels polymer glasses and thermoplastic tri- 
blocks to not explore other sets of entanglement distri- 
butions. 

Under deformation, the stress-strain relation of these 
materials is expected to follow a superposition principle, 
with one contribution from inherent glassy structures and 
their contribution to energy landscape, leading to the 
so-called yield stress, and the other contribution from 
strain-hardening effects. Events from the latter may be 
related to thermally activated processes [20l - [23| . confor- 
mational effects [28[, or topological entanglement effects 
[2^_3L^]. It was also shown in [27il'3l"| for a wide range 
of polymer glasses, that underlying entanglement effects 
contribute to the strain hardening. 

The PPA framework [H, furnishes a practical tool 
in order to identify the entanglements. Such an analysis 
has been performed on undeformed triblocks and hP and 
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Fig 2: (Color online) True stress g^z measured during bound- 
ary driven tensile test against (a) true strain A = Lz{t) / Lz{0), 
and (b) Gaussian strain g{\) — \^ ~ X^^ . Inset: Gaus- 
sian strain-hardening plot, in which stresses are normalized 
by a flow 

iP polymer glasses, averaging the number of monomers 
A^e per Kuhn segment over all chains. Results are sum- 
marized in Tab. HI in which it first appears for triblocks 
that (A^e) decreases with the increasing amount of loop 
chains, thus leading on the same manner to an increase 
of the entanglement density pe = p/2{Ne). For poly- 
mer glasses, the entanglement density pe is higher for the 
isotropic iP glass than for the directed hP glass. This 
is also expected regarding to the average imposed chains 
conformations for the latter. 

In order to study more precisely the evolution of entan- 
glements during the tensile test numerical experiment, 
the PPA has been used every 10^ MD steps on dumped 
configurations, and 10'* times during the whole deforma- 
tion process. For each dumped configuration to which 
PPA has been applied (hereinafter referred to PP chains), 
three quantities are extracted: 

[1] for each PP polymer chain, the curvature between 
a monomer i and its neighbors j = i — S/2 and 
k = i + S/2 is computed. The spatial position 
and monomer Id where the bending angle is max- 
imal is then stored. We refer to this ensemble as 
MaxBend, and 5 = 4 is used in the following. 



[2] according to a distance criterion r; = fi.rc, each 
monomer of a chain 1 having another monomer of 
a chain 2 at a distance less than ri is an element of 
the ensemble Contacts for chains 1 and 2. 

[3] a more restrictive ensemble per chain is built as 
the intersection of both previous ensembles. En- 
tanglements taking part of this ensemble can be 
considered as Mutual entanglements. 

Both Contacts and Mutual ensembles depend on the 
value of the contact length r; = p^.Tc, with p, £]0,oo[ 
and Tc the cut-off distance in Eq.([T]). Several values were 
tested. For ii < r~^, the distance criterion to build con- 
tacts is too restrictive, and the ensemble contains too few 
elements to perform a correct statistical analysis. In the 
limit of /i — >■ oo, one has "Pa+Za- (m? ^o) ~^ Oi'^{t,to) 
(see Sec. IIV A ip . as expected in this mean- field limit. 
We found optimal results ioi p = r^^ , and checked that 
varying it do not alter main results. This value will be 
used in the following. 

Finally, from a computational point of view, the com- 
plete analysis, including parallel implementation of PPA 
on a 16 processors computer and well optimized post- 
processing codes, took around 1300 s CPU per dumped 
configuration. 

III. RESPONSE TO DEFORMATION 

A uniaxial tensile test is performed at a temperature 
fc^T = 0.3w^^ under anisotropic barostatic conditions 
Pxx — Pyy — 0, with an applied upper grip velocity 
woGz = lO^^aaij/TLjex- The true stress azz is moni- 
tored as a function of the elongation A = Lz{t)/Lz{0)- 
Results are shown in Fig. ^a.) for all samples and versus 
the true strain ln{X). One can separate the stress-strain 
curve in three regions: (i) an elastic regime with possible 
localized irreversible monomer rearrangements [l9|, (ii) 
a maximum ay which marks the onset of plasticity, di- 
rectly followed by a strain softening regime, (Hi) a quasi- 
ideal plastic flow a flow and a strain hardening regime at 
large deformations. For triblocks, the first stage extends 
to strain less than 5 — 7%, while for hP and iP polymer 
glasses, this one covers the region with strain less than 
2%. From this elastic regime, the Young's modulus Ey 
is extracted, with higher values for polymer glasses com- 
pared to ones for triblocks. Results are given in Tab. HI 

The ratio of Young's moduli for homopolymers to tri- 
blocks is in quantitative agreement with one observed 
experimentally 0. We also observe that Young's mod- 
ulus for directed polymer glass hP is naturally higher 
than the one for the isotropic glass iP, for a deformation 
imposed along the hP preferential chains average orienta- 
tion. Depending the percentage of loop chains, different 
stiffness are observed for triblocks in the elastic regime. 
It appears that the 100% loop seems to be weaker, while 
it is hard to conclude between other triblocks. Larger 
systems with several lamellae are needed to go further. 

In Fig. Hljb), the true stress is also plotted in the 
Gaussian strain-hardening framework using the relation 

= o-fiow + Crg{X), where g{X) = A^ — A~*. Fitting 
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g(^) 

Fig 3: (Color online) Mean second Legendre polynomial P2 
between bond segments and applied deformation director 
upon deformation. (P2) is plotted for triblocks, for chains in 
the rubbery phase, and on (inset), for chains in the glassy 
phase. Both plots are against the Gaussian strain g{X) — 
\^ — X~^. Also plotted is the Legendre polynomial for hP and 
iP homopolymers. The increase of (P2) mimics the trend of 
polymer chain segments to align along to e,, the direction of 
applied deformation. 



the stress-strain curves with this expression defines the 
hardening modulus Gr and the offset yield stress a flow 
Results are summarized on Tab. U for all systems. This 
kind of plot allows us to separate the plastic flow and 
stress-hardening regions for 5(A) > 1, to the preyield one 
for g{X) < 0.8. By this way, it also appears for triblocks 
that stress-strain curves seem to show two peaks around 
the yield zone. 

In Figs.UJa) and (b), we note different behaviors under 
deformation for triblocks and pure polymer glasses. The 
directed hP glass displays a stronger strain-hardening 
than the isotropic iP glass. Then polymer glasses have a 
stronger strain-hardening than triblocks. For triblocks, 
large stress drops appear suddenly and earlier when the 
amount of loop chains increases. 

The presence of the elastomer also affects the response 
of triblocks. Instead of breaking, drops are nearly fol- 
lowed by local strain hardening zones, until another large 
drop appears, thus relaxing the stress in the whole mate- 
rial. This has an impact on chains orientation, mainly in 
the glassy phase, where a chevron- type structural change 
may occur 

The orientation process of chains upon deformation 
can be monitored by the second Legendre polynomial 
P2 = ^ (3cos^(0) — 1), where 9 is formed by the vector 
between monomers + 1] of the same chain and the 
direction of the applied deformation. The average orien- 
tation of the segments is taken as the ensemble average 
(P2), where (P2) = 1 means a tendency of chains to ori- 
ent along the direction of the applied deformation, and 
(P2) = means an average random orientation of chains. 

In Fig. [31 (P2) is plotted, for triblocks, for monomers 
in the rubbery phase. Average chain orientation becomes 



weaker with the increase of the amount of loop chains. 
This is expected because of the average isotropic confor- 
mations of these chains. In the inset of the Fig. |3l an 
orientation process is activated in the glassy phase, un- 
til the yield point is reached around g{X) ^ 1. After 
yield, glassy chains adopt an average different orienta- 
tion that does not change anymore and with a strength 
that increases with the amount of loop chains. This tilt- 
ing process is associated to a topological transition of 
conformations of glassy chains, due to the presence of 
discontinuities in the elastic medium of triblocks at the 
interfaces. The energy transferred to the material under 
the applied deformation leads to a localization of stresses 
(and/or specific monomer displacements) in the vicinity 
of the interfaces. This energy transferred to the glassy 
phase is strong enough after yield stress to cause an in- 
stability of the lamellae, leading to a transition from a 
lamellar to a chevron type structure . 
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Fig 4: (Color online) Spatially resolved distributions of entan- 
glements from each ensemble (see Sec. Ill Cl for details) during 
deformation. Black lines represent the average position of the 
interfaces S/B (for z < 0) and -8/5* (for z > 0), and symbols 
are related to the considered ensemble, (a) For 0% loop and 
(b) for 60% loop triblocks. 

We note that the 30% and 60% loop triblocks have a 
tendency to orient more in the glassy phase than the 0% 
triblock. Thus the tilting process is more pronounced 
when the amount of loop chains increases, meaning that 



6 



more energy is transferred trough the interface to the 
glassy phase. For the special case of the 100% loop, large 
fluctuations are observed due to the higher capacity of 
chains to reorient. Furthermore, it clearly appears that 
the changes in trend of (P2) in the rubbery phase can 
be well mapped onto the large stress drops in the stress- 
strain curve of Fig. ^h). Hence, these ones find their 
origin in the release of topological constraints and average 
chain orientation. 

In Figs, m^a) and (b), the time and space analysis 
of entanglements from ensembles Max Bend, Contacts 
and Mutual, are plotted during the whole deformation 
process. In these figures, the black lines represent the 
average position of the interfaces between S/B phases 
for z < 0.0, and B/S phases for z > 0.0. For the 0% 
loop triblock in Fig. HJa), MaxBend entanglements are 
mostly localized in the glassy regions, and in the vicinity 
of the interfaces. For the 60% loop triblock in Fig. UJ^b), 
another contribution appears in the rubbery phase, un- 
til the stress drop located at g{\) 2 (sec Fig. ^ is 
reached and after which MaxBend entanglements local- 
ize on both parts of the interfaces. Contact entangle- 
ments probe the dynamics of relevant contacts between 
chains upon deformation. For both triblocks, relevant 
contacts disappear in the vicinity of the upper moving in- 
terface, nearby to the value of the yield point at g{X) ~ 1. 
This lost of the interchain connectivity affects the balance 
with maximum curvature entanglements, then leading to 
a local interface morphology change, i.e. the chevron- 
type orientation transition. 

For the 0% loop triblock, these contacts persist in the 
rubbery phase and far from the moving interface. On 
the other way, they tend to be spatially distributed at 
a fixed distance from the moving interface for the 60% 
loop triblock, and coincide with some maximum bend- 
ing entanglements. This might be an effect of the loop 
chains, which preserve the material stiffness and stress 
transmission by adapting their conformations. It ex- 
plains why strain-hardening for 30% and 60% loop tri- 
blocks is weakly dependent of the amount of loop chains, 
in regard to the 0% loop triblock. On the other way, it 
is found (not shown) for hP and iP homopolymers that 
the same entanglement sets are well distributed over the 
whole simulation box. 

Finally, we note that the ensemble Mutual, which is 
more restrictive, seems to be a good compromise in order 
to capture relevant entanglements. In the following, wc 
will mainly focus on that ensemble. 



IV. ENTANGLEMENTS STATISTICS 

Using statistical ensembles for entanglements defined 
in the previous Sec. Ill C[ a specific attention is paid 
on the "participative dynamics" of these entanglements 
upon the tensile test numerical experiment depicted in 
Sec. mil To this aim, we split the contribution of these 
entanglements in two main contributions: from the initial 
distribution, and from the distributions on the fly. 




Fig 5: (Color online) Relative ratio (/i, t, to) of created 

vs annihilated entanglements during the whole deformation 
process. Events are taken from the initial Mutual entangle- 
ments set. Main panel: results for triblocks. Inset: same 
ratio for the iP and hP polymer glasses. The dashed line 
represents the yield point for these glasses. 

A. Contribution of initial entanglements 

The initial distributions of entanglements form the ini- 
tial networks of triblocks and polymer glasses. For their 
representation, we define the matrix Jl^ of distributions 
of events from the Mutual ensemble set at time t. This 
matrix has M rows and each row represents a chain of 

monomers. Element of a row is equal to when the 
monomer i is not an clement of the considered ensemble 
and 1 otherwise. 



1. Creation vs annihilation 

We first define a~^{iJ,,t,to) as the number of events 
present at time t but not at time to divided by the size 
of the ensemble set at time t times the amount of time 
(t — to), thus characterizing the rate of created events at 
this time. Using the above representation: 

^"^^'*'^°) = (d^^Ka-"aJ (3) 

where Nt is the size of the ensemble set at time t and 

contains complementary elements (Boolean NOT) 

with respect to the universe (0,1) for each state defining 

a monomer in the transposed matrix '^ft^ of . Hence, 
^ _ — — ^ 

the operation O : £t — > 11 is related to the Boolean 
NOT transformation applied to elements of each row. In 
practice, if a monomer of a chain has a state in a row 
of rj at to, then such a monomer will have a value 1 in 

— *o ^ 

the column of the associated transposed matrix ^ft . 

— to 

Conversely, we define a~ {fj,,t,to) as the number of 
events present at time to but not at time t and also di- 
vided by the size of the ensemble set at time t times the 
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time length (t — to), thus characterizing the rate of an- 
nihilated events at measuring time. In a more compact 
notation, one has: 

a-(u,i,to) = - 3_ — Tr(h ■'^U ] (4) 

' {t-to)Nt V=t =toJ ^ ' 

where Jl^ contains complementary elements (Boolean 
NOT) with respect to (0, 1) for each state defining a 
monomer in the matrix 17^. Thus, if a monomer of a 
chain has a state in a row of 11^, it will become 1 in the 

same row of £2^ . 

The quantities Eq.([3|) and Eq.(|3]) may differ during 
the deformation, leading to an asymmetric creation- 
annihilation process. The ratio "Pa+Za- (Mj ^o) captures 
such a process. 

In Fig. [SJ the variation of the ratio Va+ /a- (/^i ti to) is 
plotted for the Mutual ensemble set of entanglements. 
For triblocks, and before yield at g{X) ^ 1, balance be- 
tween alive and created events is achieved with regard 
to the initial entanglements set. Then after yield, more 
and more new events are created, which unbalances the 
ratio with annihilated ones, then leading to an asym- 
metric creation-annihilation process of entanglements. It 
coincides with the chevron-type transition in the inset 
of Fig. [3] that allows glassy chain ends to release some 
topological constraints. The asymmetric process has an 
higher magnitude in the case of the 100% loop sample, 
with changes in slope that also correspond to the large 
drops appearing in the stress-strain curve in Fig. [2] 

After yield, these more created or activated events can 
be seen as the result of "cascades" of disentanglements 
of older events, where such processes alter chain confor- 
mation and, this, in order to maintain the rigidity of the 
system. It means that the 100% loop triblock is more af- 
fected by the deformation in the plastic flow regime than 
other triblocks and needs to explore more chain confor- 
mations in order to preserve a mechanical stability. 

Intriguingly, it appears that the 30% loop triblock fol- 
lows the same trend when compared to the 60% loop one. 
The latter is indeed closer to the 0% loop behaviour. A 
lower balance is then achieved for this triblock, which 
means that it is not forced to explore more and more 
new sets and destroy older ones in order to preserve its 
mechanical stability. In that sense, this amount of loop 
chains seems to much preserve the triblock integrity, and 
plays in favour of the future metastable states during the 
deformation. This can be related to the more pronounced 
chevron-type process for this triblock, which compen- 
sates the loss of topological constraints due to the low 
amount of crossing chains. A this point, a more precise 
study is needed to explore this compensation effect, and 
how this one depends on the amount of loop chains. 

In the inset of Fig. [SI the ratio T'a+ /a- (/^i h) is plot- 
ted for the iP and hP polymer glasses. Before yield, 
balance between created and annihilated events is ob- 
served, followed after yield by an asymmetric process. 
Its magnitude for the iP glass is higher than the one for 
the directed hP glass, because new metastable topolog- 
ical intermediate configurations are more accessible for 
iP. Despite of the glassy state, and under deformation, 




Fig 6: (Color online) Distribution 'Ps(Mj*i^o) of survival 
events between initial undeformed state to and deformation 
time t. Events are part of the Mutual ensemble set of entan- 
glements, and data are plotted for triblocks. The 30% loop 
chains triblock is omitted for clarity. Inset: same quantity for 
hP and iP polymer glasses. The vertical dashed line marks 
the onset of plastic flow. 

energetic barriers have then an higher probability to be 
crossed, leading to this higher ratio. Finally, the ratio 
Pa+ /a- (Mj to) for iP and hP glasses is also weaker than 
the one for triblocks. It is mainly due to the presence of 
a rubbery region in triblocks, in which chains can explore 
more conformations. 



2. Survival events 

Additionally, we define the distribution VsifJ-Tt^to) of 
Mutual survival events between initial time to and the 
measuring time t during deformation. Using the same 
formalism for Eqs.dS]) and (|3]), it can be written: 

Vs{f^,t,to) = ^Tr[n-^n^) (5) 

Such a distribution is plotted in Fig. [51 It clearly ap- 
pears two distinct regimes with increasing g{X): a plateau 
until yield strain is reached, followed by a power law de- 
crease of survival events Vs{lL,t,to) ^ g{X)~^ . For tri- 
blocks, 13 « 0.83, 0.74, 0.68, 0.97 for 0, 30, 60 and 100% 
loop chains respectively. Hence, for the three first tri- 
blocks, the plateau value before yield is quantitatively 
the same. Except the 100% loop triblock, the rate of 
loss of the contribution of initial mutual entanglements 
after yield decreases with the increase of the amount of 
loop chains. The triblock with the lowest amount ex- 
plores more rapidly new entanglements set on deforma- 
tion than ones with a higher amount; loop chains offer a 
longer memory effect of the initial mutual entanglements 
network. 

The plateau can be interpreted as the fact that ini- 
tial Mutual entanglements mainly contribute to the elas- 
tic linear (and nonlinear) part of the strain-stress curve. 
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while their contribution to the plastic flow and strain- 
hardening regime decreases with increasing strain. An 
higher preyield "plateau" is observed for the 100% loop 
chain triblock, in which chains more preserve the ini- 
tial network of mutual entanglements, but on counter- 
part make the material stiffer and more sensible to high 
deformations. An higher postyield power-law exponent 
is also measured for this triblock, which means that this 
material needs to explore new states more rapidly. 

Finally, in the inset of Fig. [6j Eq. ([5]) is plotted for 
the hP and iP glasses. Due to their chain average ori- 
entation, the preyield regime depicts a higher survival 
magnitude for the directed glass than for the isotropic 
one. In the postyield regime, the power law behavior is 
fitted with an exponent j3 « 1.01 for hP polymer glass, 
and /3 « 0.2 for iP glass. The rate of loss of initial mu- 
tual entanglements is then higher for the directed glass 
than for the isotropic one. It is because, given the ini- 
tial hP glass entanglement network, the chain average 
orientations are too constraining to expect a mechani- 
cal equilibrium in the postyield regime. This glass needs 
to explore new mutual entanglement states more rapidly 
upon deformation than the isotropic one. 
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Fig 7: (Color online) Distribution Vsifi, A*) of survival events 
between deformation times t and t—A*. Events are part of the 
Mutual ensemble set of entanglements, and distributions are 
plotted for triblocks with specified percentage of loop chains. 
The parameter A* is given in tlj units. Exponents a and 
7 that characterize the two regimes discussed in the text are 
also drawn. The vertical dashed line marks the yield point. 



B. Contribution of on the fly entanglements 

The participative dynamics of mutual entanglements is 
evaluated by monitoring the contribution of their distri- 
butions to the current deformation time t, given a span- 
ning time A*. To this aim, we define the retarded distri- 
bution VsilJ-, A*) of survival events between time t and a 
new origin to ~ t — A*. This distribution can be written 
using the same formalism as for Eq.©: 

P,(A.,A*) = i-Tr(£^-^£^,) (6) 

In Eq.®, the quantity is then related to the 

matrix of distribution of mutual entanglements at a given 
origin A*. The distribution Vs{n, A*) is plotted in Fig. [7] 
for the different SBS triblocks and their amount of loop 
chains, i.e. from 0% to 100%, versus the Gaussian strain 
g{X) = - A-i and A = L,{t)/L,{0). 

In Fig. [71 the "probe" times A* are given in tu 
units, with all A* < ttot/2 and the total amount of 
simulation time ttot = 6 x lO^r^j for all simulated sys- 
tems. Values of the probe times are taken according to 
specific deformation zones of the stress-strain curves in 
Fig. [2j For the elastic zone, where teiasUc is an estimate 
of the upper time for the elastic limit, we have A* = 
1.2 X 10^, 1.2 X 10'' < teiasuc- The value A* = 6 x 10'' 
is chosen to lie between tpi and tp2, which are estimates 
of the deformation time associated with the two peaks in 
the yielding zone of the stress-strain curve. Finally, both 
values A* = 1.2 x 10^, 1.8 x 10^ > thard, where thard is 
the time after which strain-hardening occurs. 

It appears in Fig.[7]before yield and A* < teiasUc, that 
distributions have a lower magnitude when the amount of 
loop chains increases, except for 100% loop. If loop chains 
maintain the initial set of entanglements, the 0% triblock 



is more stable than other triblocks with regard to new on 
the fly entanglement sets. The 100% loop triblock depicts 
higher distributions, which is a direct effect of the special 
amount of these chains as already discussed in Sec. lIV Al 

The 'Ps(m, A*) curves seem to follow a master curve 
which can be split in three parts: (1) a plateau for A* < 
^elastic, which stays until the yield strain at g(A) 1 
is reached, (2) a. power law decrease fitted with an ex- 
ponent a, and (3) a, power low increase also expressed 
in terms of an exponent 7. The transition from the a 
to 7 regime gives rise to a minimum, which we denote 
min{7',(M,A*)}. 

Events in the a regime are related with initial events 
from the preyield zone, and quantify their contributions 
to events moving away from the yield point. The min- 
imum 77im{7's(/^, A*)} is related to the contribution of 
events at yield point to postyield ones. Finally, the 7 
regime characterizes the contribution of events from the 
post-yield region to ones in the strain-hardening regime. 
Hence, an increase of the contributions in this regime 
marks the emergence of long-time correlations and mem- 
ory effects. 

In the 7 regime, the distributions of survival mutual en- 
tanglements separated by the shortest A* have a higher 
magnitude than ones for longer A*. On counterpart, the 
route to the memory effect is less pronounced for small 
probe times A* than in the case of larger probe times if 
we relate this route to the 7 exponent. Thus, at larger 
probe times, survival distributions between postyield and 
in- strain-hardening entanglements tend to higher magni- 
tudes more "rapidly" upon deformation. Finally, these 
trends in the 7 regime are also more and more pro- 
nounced when the amount of loop chains increases, as it 
is shown in Fig.|SlJb), where the "memory" exponent 7 is 
plotted against the probe time. Consequently, memory 
effects and long time correlations between sets of mu- 
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Fig 8: (Color online) Variation of a (a) and 7 (b) exponents 
depicted in Fig. [7] against the probe time A* that separates 
two measuring times upon deformation. Exponents charac- 
terize the two power-law behaviors for Vsi^i, A*). On both 
main panels, data are plotted for triblocks, and on insets for 
hP and iP polymer glasses. 



tual entanglements are more and more important in the 
postyield regime, and the loop chains reinforce the mem- 
ory effects. 

In Fig. Ela), the increase of \a\ with A* denotes an 
increase of the loss of the contribution of preyield mu- 
tual entanglements to postyield ones. Thus, long-time 
correlations become weaker and weaker between the sets 
of mutual entanglements from these two regimes. The 
values of |a| for the 100% loop triblock are also higher, 
especially for A* < tpi.p2, which points out that this tri- 
block experiences more chain conformations than others. 

Exponents a and 7 for the polymer glasses are plotted 
in the insets of respectively Figs, ^a) and (b). For the 
isotropic glass, the loss exponent a is higher than for the 
directed glass. It is because intermediate conformation 
states are more accessible for the isotropic glass, and thus 
it can explore more mutual entanglement sets. The slope 
of |a| with respect to A* is also lower when compared 
to triblocks, because of the pure glassy state that slows 
down the debonding. 

Finally, the memory exponent 7, which characterizes 
the postyield to strain-hardening contributions of sets of 



mutual entanglements, depicts, in Fig. [S^b), a common 
trend for both polymer glasses. With the increase of 
the probe time A*, it then appears that memory effect 
and long-time correlations become stronger. This trend 
seems to be more pronounced than for triblocks, which 
is again due to the pure glassy state. 



V. CONCLUDING REMARKS AND 
OUTLOOKS 

During this study, an isotropic, a directed poly- 
mer glass and lamellar directed triblock copolymers 
were generated. The directed triblocks were devel- 
oped in order to mimic the mechanical behavior of 
glassy /rubbery/glassy nanostructurated materials in the 
strong segregation regime. By tuning the amount of un- 
crossing or loop chains, it was possible to simulate dif- 
ferent glassy /rubbery interface definitions. Material me- 
chanical properties were investigated using a boundary 
driven tensile test experiment. 

Using this methodology, the impact of loop chains 
on the mechanical response of triblocks has been ap- 
proached. Their stress-strain curves depict large stress 
drops that come earlier upon deformation, with the in- 
crease of the amount of loop chains. A chevron- type tran- 
sition is also observed, with a strength that increases 
with the same amount. Besides this facts, the global 
mechanical properties seem to not be deeply altered by 
these amounts, so that, it is hard to predict a critical 
amount of loop chains above which triblocks integrity is 
achieved. To overcome this fact, larger systems with re- 
peated lamellar units should be considered in order to 
improve statistics, decrease fluctuations and reach a self- 
averaging regime. 

The sets of mutual entanglements between chains were 
monitored during deformation, using a primitive path al- 
gorithm. A specific attention has been paid to the study 
of their statistics. 

Given the initial undeformed sets of mutual entangle- 
ments, the evolution of the ratio of created over anni- 
hilated events, as well as the distributions of survival 
events, has been monitored. It has been shown that 
these quantities follow the two main parts of the stress- 
strain curve, namely, the preyield and postyield regimes. 
Monitoring on the fly sets of entanglements during the 
deformation allowed us to also approach the internal 
metastable cross-link states. 

It appears that during the strain-hardening regime, the 
debonding of the initial mutual entanglement network 
occurs. Debonding of on the fly entanglement sets is 
correlated in time, depending the initial time window 
from which depart these sets. If one considers sets from 
the preyield regime, their contributions to the postyield 
regime decrease with a power-law behavior denoted by a 
loss exponent a. On the other way, ones that depart from 
the postyield regime depict memory effects in the strain- 
hardening regime. In that case, a power-law behavior is 
observed, which is characterized by a memory exponent 
7. The same trends occur for polymer glasses as well as 
for triblocks. 
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For triblocks with a non-null amount of crossing 
chains, the lower the amount is, the longer the memory 
effect of initial undeformed sets of mutual entanglements 
is. On the same way, loop chains also enhance long-time 
correlations and memory effects of internal metastable 
cross-link states continuously explored upon deforma- 
tion. 

Finally, it should be noted that the results obtained 
may depend on the applied strain rate. Further works 
should be employed in order to quantify this dependence. 

The methodology developed during this study gives 
rise to further insights. First, the created and activated 
sets of mutual entanglements emerge from a complex pro- 
cess, which could contribute to new stable intermediate 
conformations as well as to participate to destabiliza- 
tion processes and related unstable states. The inner 
mechanisms of propagation of these ensemble sets should 
be more precisely studied, and especially how the large 
stress drops in Fig. [2] could be interpreted in terms of 
entanglements and how the amount of loop chains com- 
pensates the loss of topological constraints, as observed in 
FiglS] With the methodology developed in this study, one 
could access these localized or delocalized (long-range) 
processes. A quasistatic protocol could be suitable in 



order to remove the thermally activated processes and 
distinguish the transmission mechanisms. 

Furthermore, this methodology allowed us to deal with 
the problematic of the way external forcing is "inter- 
nalized" by the systems via the cross-link process. The 
rheological response of these systems involves permanent 
cross-linked and uncross-linked sets, which can be de- 
scribed by some surviving probabilities C(^)j until t times 
of deformation occur after cross-linking. These states 
hold internal stresses. Under external forcing, they are 
metastable and require far from equilibrium operations in 
order to be created. In the studied systems, which depict 
rubbery phase or "interchain" connectivities, the gener- 
ation of internal states as metastable states is closely re- 
lated to the notion of plasticity. This notion is itself re- 
lated, probably among many others, to some time scale 
competitions that take place during this generation pro- 
cess. This may lead to an instantaneous momentarily 
decrease of time scales, leading to a bias in favor of the 
future metastable state. Then, it would be also inter- 
esting to extend the notion of entanglement sets to the 
study of internal stress dissipation mechanisms and local 
displacement fields. 
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